Photon gas dynamics in the early universe 



Chris Clarkson^'^'B and Mattias Marklund^'0 

^ Cosmology and Gravity Group, Department of Mathematics and Applied Mathematics, 
University of Cape Town, Rondehosch 7701, Cape Town, South Africa 
'"^Institute of Cosmology and Gravitation, Mercantile house. 
University of Portsmouth, Portsmouth, POl 2EG, United Kingdom 
'^Department of Physics, Umea University, SE-901 87 Umea, Sweden 
(Dated: February 2, 2008) 

Quantum electrodynamics predicts that photons undergo one-loop scattering. The combined effect 
of this on the behaviour of a photon gas for temperatures above ~ 10^° K results in a softening of 
the equation of state. We calculate the effect this has on the effective equation of state in the early 
universe, taking into account all the species of the Standard Model. The change to the dynamics of 
the early universe is discussed. 
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Introduction The quantum vacuum has been a topic 
of fundamental interest in physics since the advent of rel- 
ativistic quantum theory, and has spawned a number of 
experimentally verified concepts, of which a few are the 
Lamb shift, Delbriick scattering, photon splitting, and 
pair creation [Ij]. In particular, the possibility of probing 
its structure using state-of-the-art experimental facilities 
has recently been the in the focus of a number of works 
(see e.g. [2| and discussions therein). It is well-known 
that the properties of the quantum vacuum can have im- 
portant consequences for black hole (BH) physics Q , in- 
troducing the concepts of BH thermodynamics [3| . The 
properties of the quantum vacuum also affects the prop- 
agation of photons, introducing a background dependent 
phase velocity v of pencils of light, such that f ^ c. This 
is a result of so called elastic photon-photon scattering 
in which virtual electron-positron pairs mediate the col- 
lisions between photons. Most of the discussions in the 
literature regarding the effects of background fields on 
the null-cone structure has been limited to the weak field 
case. This is relevant when possible experimental setups 
are discussed. However, the propagation of photons on a 
very strong background field yield interesting new effects. 

We shall consider the full one-loop quantum electro- 
dynamic (QED) correction [given by the lagrangian ([1]) 
below] represented in terms of Feynman diagrams as 
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In Ref. '6] it was shown that this contribution leads to a 
nontrivial implicit dispersion relation for intense incoher- 
ent monochromatic photons. For super-critical intensi- 
ties, i.e., for radiation energy densities > eojii^critP, where 
Ec.it = m^C'leh k. 1.32 x lO^^V/m {h, is the Planck 
constant, rrie is the electron mass, and c is the vacuum 
speed of light) , is the critical field, this leads to a notice- 
able decline in the photon phase velocity which for very 
high energy densities eol^^critP) approaches the value 
V ~ cjyfh. Also, since the phase velocity is reduced, col- 



lective photon pair creation through the Schwinger mech- 
anism is not possible even though the gas intensity ex- 
ceeds the corresponding Schwinger limit. It should be 
noted that for a very hot thermal photon gas, a signifi- 
cant part of the Planckian distribution will be above the 
Conipton frequency meC^/?i, thus making single photon 
pair production viable. In a cosmological context, the re- 
sult of such a nontrivial dispersion relation is to give an 
effective equation of state (EOS) of radiation Q of the 
form -p^ — p-y/Sn^, where p is the isotropic pressure, 
is the radiation energy density, and n — cjv ^ 1 is the 
refractive index of the radiation (see also |8| ) . 

In this paper we shall consider the weakening of the 
photon equation of state on the evolution of the early 
universe. At the onset of primordial nucleosynthesis, the 
effect of the change in refractive index on the equation 
of state is small, but as will be shown below such reduc- 
tions in the effective radiation pressure at earlier stages 
introduce noticeable effects. 

The photon gas equation of state The vacuum elec- 
tromagnetic lagrangian for arbitrary field strengths at 
the one-loop level is given by the Heisenberg-Euler la- 
grangian 0] 
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which is valid when the photon frequency satisfies lo <C 
c^e = m^c^/n. Here, a = rc-^s ^2^1/2 , . 1/2 
1 1/2 



where ^ = \FabF'''' = 
'^"^FnhF,^ = -cE-B, a is the fine- 



2(c2B2-E2),^^ = ie--w,d 
structure constant, Cq is the vacuum permittivity. The 
critical energy density is pcrit = eol-E-critP ~ 1.55x10^^ kg 
m~^s~2 = (7r2fc^/15?i^c'^)T^^j.ij corresponding to a critical 
temperature Tcdt « 1-2 x 10^°K « 1.0 MeV. Hereafter, 
we shall use units where Ti — ^ c = I- 

As shown in the equation governing the refractive 
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index as a function of = p-y/ Pcrit is given by 
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where the function F is 
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where the integral is taken on a contour to the right of 
the imaginary axis. The solution to Eq. Q may be 
obtained by considering the low energy approximation, 
F{x) ~ il^;^, which will be sufficiently accurate for the 
effects considered here. As noted below Eq. (1), the gov- 
erning Lagrangian is valid for arbitrary field strengths, 
but with the photon frequency lo uje- For a ther- 
mal distribution, a significant number of the photons will 
have a frequency uj ^ uje when T ~ Tdit. Thus, the 
Planckian photon distribution implies that the present 
model i formally valid for T < Tent- However, in what 
follows we will still make use of the Lagrangian (1) for 
temperatures > Tdit in order to probe the consequences 
of such elastic photon-photon scattering. The high-T 
model for this process, in a cosmological context, is a 
subject for future studies. We note that the refractive 
index satisfies 1 < n < y/E, as may be seen from Eq. ([2]). 

Using the low energy approximation, the equation of 
state parameter for a photon gas, u'-y(C) = PiIp^ = 
l/3n^, is the solution of 
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Since C, = (T/TcHt)^, we see from the small dimensionless 
constant on the right hand side that photon-photon scat- 
tering will start significantly affecting the dynamics of a 
photon gas around T ~ lOTcrit ~ lOMeV. However, at 
such high temperatures, a black body spectrum will give 
rise to a significant pair production. The rate of such 
pair and photon production can roughly be modelled ac- 
cording to c and 

tively. Here aj^-ep and (Tep_^^ are the cross section for 
pair production and pair annihilation (in the center-of- 
mass frame), respectively, and and Ue are the photon 
and electron number densities, respectively. Such direct 
pair production by high frequency photons will imme- 
diately result in a nonlinear cascading process, dissipat- 
ing the electromagnetic energy into fermionic degrees of 
freedom. It is possible to model this process more accu- 
rately using kinetic models, taking higher order interac- 
tions 0. However, a simple, but descriptive, local model 
is given by dtn^ « i/^n^ — Vepn,,, dtn-e ~ —v-yn-y + VepUe, 
that gives rise to a typical initial exponential growth in 
the electron number density which reaches a nonlinear 
equilibrium at later times. Furthermore, the possibil- 
ity of multi-photon pair production could be accounted 



for in a more accurate model. However, the time-scale 
for pair-production for a high-temperature 1 MeV) 
thermal photon distribution is short on a cosmological 
time-scale, and for our purposes it suffices to introduce a 
cut-off for the photon energy at 1 — 10 MeV (see Ref. Q 
for a discussion of the thermal range for pair production). 
Moreover, the onset of thermal equilibrium between pho- 
tons and electron-positron pairs has only recently been 
investigated in detail using kinetic theory For small 
temperatures, however, we have 
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which is actually a good approximation to ~ 5MeV. 

Contribution to the overall EOS In the early universe, 
at temperatures above roughly one MeV, photons are no 
longer a major contributor to the overall energy den- 
sity. As the temperature is increased above the rest 
mass of particles in the Standard Model they become 
relativistic, contributing to the energy density and pres- 
sure through the effective number of degrees of freedom, 
g*{T) = (30/7r2)(pt,ot/r^), where ptot is the total energy 
density, and T is the temperature of the photons [l3| . 

The total energy density at temperature T is given by 
a sum over all species in equilibrium in terms of phase 
space integrals [10|, 
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while the total pressure is given by 
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Here, the — sign refers to bosons, and the +-sign to 
fermions. We have used xi — mi/T^ and gi is the num- 
ber of helicity states of the i'th species. (We are as- 
suming the chemical potential is negligible.) Normally 
in the expression for the pressure the photons may be 
included in the sum; here we have written it separately 
in order to use the modified EOS, Eq. The energy 
density of the photons as a function of temperature may 
be calculated from the appropriate term in Eq. ([6]), i.e., 
Py = [tt^ /lb)T'^, while the modified pressure is given by 
the effective QED equation of state — Wyp-y. 

For some species the temperature Ti may not be the 
same as the photons, while still having a thermal distri- 
bution. Here, we consider only the effects of neutrinos 
having this property: for temperatures above a few MeV 
the neutrinos and photons share the same temperature. 
When the electrons and positrons freeze-out their entropy 
is transferred to the photons, raising their temperature 
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relative to the neutrinos by = (4/ll)y^T, which may 
be found by conservation of entropy [lO|.^ 

We now define the effective equation of state parameter 
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which is of course equal to 1/3 for purely relativistic 
species. However, as particles freeze out and become non- 
relativistic, the contribution they make to WeS causes it 
to deviate from 1/3 during the freeze out. We shall anal- 
yse the effective equation of state as a means of quanti- 
fying the QED correction to the photon gas equation of 
state. We do this by first considering WeS in the standard 
case of ui^ = 1 /3 and then compare it to as given by 
Eq. Q. 

The difference between the equations of state is then 
given by 
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For low temperatures we find the correction to the stan- 
dard EOS as 
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For temperatures around 1 MeV we have 5, « 10.6, and 
therefore find that the EOS is reduced by around 0.01% 
at r = l.SMeV due to photon-photon scattering. This 
reaches 1% by T ~ 5MeV, where we must take into ac- 
count the non-linear equation of state for given by 
Eq. (UJ; by 10 MeV, when the photons are contributing 

18% of the total energy density we find that the EOS 
is reduced by more than 10%. Of course, it is unlikely 
that the model we have used is sufhciently complex at 
such high temperatures. 

Dynamics of the early universe Given our effective 
equation of state as a function of temperature we may 
now calculate the dynamics of the early universe. Be- 
cause the density as a function of temperature is un- 
changed by photon-photon scattering, the Hubble rate as 
a function of temperature will be unchanged too, since 

The behaviour of the scale factor a, and time, t, with 
temperature will be different between the two cases, since 
Ptot — ~3i?/9tot(l + Wcff). Obviously this must be calcu- 
lated numerically in general. To demonstrate roughly 
how photon-photon scattering modifies the evolution, let 



us define e(T) = 1/3 — WeS > 0. It may be shown that 
e ^ 0.05 whatever EOS we use for the photons. To first- 
order in £ we find that the scale factor behaves as 
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The time evolution with temperature, on the other hand, 
is 
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where k — 8ttG. Therefore, we may estimate the low 
temperature change to the time and scale factor evolution 
to leading order as 
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where we have assumed that is constant. The evolu- 
tion of the scale factor evolves as normal at late times 
(low temperature). Given the time difference between 
the two cases, the time between, e.g., the quark-hadron 
phase transition and the onset of nucleosynthesis will be 
shorter when photon-photon scattering is taken into ac- 
count. 

Discussion We have considered the collective effect 
that photon-photon collisions (currently under experi- 
mental probing [2]) have on the dynamics of the early 
universe, with surprising results. Taking into account the 
full one-loop Lagrangian for photon-photon scattering on 
an incoherent radiation gas results in a drop in pressure 
of the photons. While the change in EOS is small at 
the onset of nucleosynthesis, the change in the dynam- 
ics of the universe before this time could have an effect 
on high-energy phenomena such as primordial black hole 
production 13]. However, more precise estimates of such 
effects will require further research. 

Photon-photon scattering should therefore be built 
into models of the early universe in the hadron era and 
later. Indeed, as primordial nucleosynthesis calculations 
are now accurate to 0.1% [l^, the fact that the neutrino 
temperature is increased by this fraction imply that a 
further refinement may be necessary in these models. 
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^ It may be shown that the neutrino freeze-out temperature is 
increased by around 0.1% by the effects considered here. 
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